%test lda
clear;

c1=[1 2;2 3;3 3;4 5;5 5];
c2=[1 0;2 1;3 1;3 2;5 3;6 5];


%using LDA
data=[c1;c2];
target=[zeros(length(c1),1); ones(length(c2),1)];
w=LDA(data, target);

u1=mean(c1);
u2=mean(c2);
s1=4*cov(c1);
s2=5*cov(c2);
sw=s1+s2;
sb=(u1-u2)'*(u1-u2);

v=inv(sw)*(u1-u2)';
[v2 d2]=eig(inv(sw)*sb);
